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Abstract. We consider FLRW cosmological models with standard Friedmann equations, but 
leaving free the equation of state. We assume that the dark energy content of the universe is 
encoded in an equation of state p = f(p), which is expressed with most generality in the form 
of a power expansion. The inclusion of this expansion in Friedmann equations allows us to 
construct a perturbative solution and to relate the coefficients of the equation of state with the 
formation of singularities of different types. 



1. Introduction 

Since the end of the 20th century there has been an increase of evidence for an accelerated 
expansion of the universe pQ. These puzzling observations do not fully fit in our present 
theoretical framework and therefore there are two trends of development attempting to ammend 
it. 

One trend is grounded on modifications of general relativity as the correct theory of gravity 
[2]. The other major trend assumes the validity of general relativity and postulates the existence 
of an exotic component in the content of the universe known as dark energy [3]. 

So far the only final state for our universe could be either infinite expansion or a final collapse, 
the big crunch. But these new scenarios open the possibility of a different final fate, since both 
violate some of the classical energy conditions [4J . New candidates for a final fate of our universe 
are the big rip [5], sudden singularities [6J, big brake [7], big freeze [SJ, inaccesible singularities 
[S], directional singularities [TO], ^-singularities [TT], in braneworld models [P2] , among others. 

In [13] we dealt with such singular scenarios from the point of view of geodesic completeness 
in modified theories of gravity. In this talk we do it following the approach of a dark energy 
fluid. More details can be found in [13] . 



2. Equations of state 

Instead of restricting to a specific model, we propose a framework which may comprise most of 
them. With this aim in mind, we assume the validity of standard Friedmann equations, 

h 2 = p , p + m{ P + p) = Q, (i) 

assuming that pressure p has a power series in density p of the matter content around a value 
p*, for which a qualitative change of behaviour is expected, 

p=p (p-p*r+pi(p-p*r + ---, (2) 

where the exponents are real and ordered, cxq < a\ < ■ ■ ■. 

Combining cosmological equations with our power expansion, we get 

p = -3Vpt (p* + kp - p*) + \ {p ~ p * )2 + ■■■ + P0 (p - P *n +pi(p- p*r + •••)• (3) 

\ 2 8 p* J 

We may integrate it for p in powers of (to — t), where to satisfies p* = p(to), if it exists, 

p(t) = p, + Pl (t - if* + p 2 (t - tf* + ■■■. 

In order to compare our results with previous ones in terms of expansions of the scale factor of 
the universe [15], 

a(t) = a(t) = c \t - t \ vo + Cl \t - t^ 1 + ■ ■ ■ , (4) 
we are to relate both series through Friedmann equation, 

2 

p = tt^-ti + —voivi - m)(t Q - 1)"~ 2 + 2 -^ m { m - m )(t - 1)^~ 2 + • • • , 
[t - t y c c 

These expansions cope with almost every cosmological model, except those with oscillatory 
behaviour as the ones in |16| . The results are consigned in table [TJ 
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3. Classification of singularities 

We follow the popular classification of singularities in [3]: 

• Big bang / crunch: zero a, divergent H, density and pressure. 

• Type I: "Big rip": divergent a. [5] 



• Type II: "Sudden": finite a, H, density, divergent H and pressure. [6] 

• Type III: "Big freeze": finite a, divergent H, density and pressure. [5] 

• Type IV: "Big brake": finite a, H, H, density, pressure, divergent higher derivatives. [7] 

In [T7] a detailed analysis of the strength of these singularities is performed, making use of 
the concepts of geodesic completeness [18] and the definitions of strong singularities due to Ellis, 
Schmidt, Tipler and Krolak [19]. Sudden and big brake singularities cannot be considered as 
such [20j, since the cosmological spacetime may be continued beyond the singular event. And 
particles travelling at the speed of light do not experience the big rip, since they need infinite 
proper time to reach it [T7]. These results may be checked in table [2] in terms of the exponents 
of the expansion of the density. Conformal diagrams for these models can be found in [21] . 



Table 2. Singularities and density expansions. 
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4. Singularities and equations of state 

Now we have all the ingredients for checking the formation of singularities in cosmological models 
with an equation of state of the form ([2]). We may integrate ([3|) in terms of power series of time 
and get the relevant exponents and coefficients. Comparing them with the results in table [21 we 
are ready to obtain the types of singularities which may appear in these models. The results 
may be found in table [3l 



Table 3. Singularities and equation of state. 
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5. Conclusions 

In this talk a detailed classification of the future behaviour of FLRW cosmologies in terms of 
the equation of state is provided. It depends at most on the exponents of the first two terms of 
a power expansion of the equation of state. The scheme provides an easy route to conclude the 
sort of singular behaviour. 
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